7 The multipole expansion 



7.1 Multipole expansion of the scalar wave equation 

Consider the emission and scattering of electromagnetic radiation. This type of 
problem involves solving the vector wave equation. The solutions of this equation 
in free space are conveniently written as an expansion in orthogonal spherical 
waves. This expansion is known as the multipole expansion. Let us examine this 
expansion in more detail. 

Before considering the vector wave equation, let us consider the somewhat 
simpler scalar wave equation. A scalar field ^(r,t) satisfying the homogeneous 
wave equation 

VV-ig = (T.1) 
can be Fourier analyzed in time 

iP(r,uj)e-''^*duj (7.2) 

-oo 



with each Fourier harmonic satisfying the Helmholtz wave equation 

(V' + fc2)V^(r,u;) = 0, (7.3) 

where = cj^/c^. We can write the Helmholtz equation in terms of spherical 
polar coordinates (r, 9, cp): 



= 0. (7.4) 



As is well known, it is possible to solve this equation via the separation of vari- 
ables: 

i^{r,uj) = Y,fim(r)Yim{e,^). (7.5) 

l,m 

Here, we restrict our attention to physical solutions which are well behaved in the 
angular variables and The spherical harmonics Yijn{9, (/?) satisfy the following 
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equations: 



= rn Yim, 



1 d . d_ 1 d'^ 

sin e de ^'"^ de ^ sin^ e d^p^ 



(7.6a) 



(7.6b) 



where I is a non-negative integer, and m is an integer which satisfies the inequality 
\m\ < I. The radial functions fim{f) satisfy 



where there is no dependence on m. With the substitution 

uiir) 



(7.7) 



.1/2 ' 



(7.8) 



Eq. (7.7) is transformed into 



\ d ^2 
dr^ r dr 



(I + 1/2)' 



ui{r) = 0. 



(7.9) 



It can be seen, by comparison with Eq. (5.39), that this is a type of Bessel's 
equation of half-integer order I -\- 1/2. Thus, we can write the solution for fim{i^) 
as 

/'".W = ^ J>+i/2{kr) + ^ Y,+iMkr), (7.10) 

where Aim and Bim are arbitrary constants. The half-integer order Bessel func- 
tions J/_|_i/2(^) and Yij^i/2{z) have analogous properties to the integer order Bessel 
functions Jm{z) and Ym{z). In particular, the J/_|_i/2(^) are well behaved in the 
limit 1^1 — >• 0, whereas the y/_|_i/2(^) are badly behaved. The asymptotic expan- 
sions (5.43) remain valid when m — / -|- 1/2. 

It is convenient to define the spherical Bessel functions ji{r) and yi{r), where 

1/2 



Jl+l/2(z), 
^+l/2(^)- 



(7.11a) 
(7.11b) 
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It is also convenient to define the spherical Hankel functions 

h^^'\z)=ji(z]±iyi{z). (7.12) 

For real z, h^'^\z) is the complex conjugate of h^^\z). It turns out that the 
spherical Bessel functions can be expressed in the closed form 



In the limit of small argument 



z' 



At 

where (2/+1)!! = (2/ + l)(2/-l)(2/-3) ■ ■ ■ 5-3-1. In the limit of large argument 

. , , sin(^ - /7r/2) ^ ^ x 

JzW ^ (7.15a) 

z 

cos(^ - /7r/2) ri X 

y/(2;) ^ ^ (7.15b) 



and 



1 z 



hf) i_. (7.16) 



The inhomogeneous Helmholtz equation is conveniently solved using the Green's 
function Gi^{r,r'), which satisfies (see Eq. (2.109)) 

(V' + fc2) G^{r,r') = -S{r - r'). (7.17) 

The solution of this equation, subject to the Sommerfeld radiation condition, 
which ensures that sources radiate waves instead of absorbing them, is written 
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(see Section 2.13) 



, i k\r—r' I 



47r \r — r' 



(7.18) 



The spherical harmonics satisfy the completeness relation 

oo I 

E E Yii,{e\^')Yim{eM = - ^')K^ose - cose'). 

1=0 m=-l 

Now the three dimensional delta function can be written 



It follows that 



oo I 



1=0 m=-l 

Let us expand the Green's function in the form 

GUr,r') = Y,9i{r,r')Yf^{e',^')Yim{e,^). 

l,m 



Substitution of this expression into Eq. (7.17) yields 



2d_ 

dr'^ r dr 



/(/ + 1) 



91 = 



5{r - r') 



(7.19) 



(7.20) 



(7.21) 



(7.22) 



(7.23) 



The appropriate boundary conditions are that gi is finite at the origin and corre- 
sponds to an outgoing wave at infinity (i.e., g oc e''^^ in the limit r — > oo). The 
solution of the above equation which satisfies these boundary conditions is 



gi(r,r')=Aji(kr<)h^'\kry), 



(7.24) 



where r<; and r> are the greater and the lesser of r and r', respectively. The 
correct discontinuity in slope at r = r' is assured if A = ik, since 



(7.25) 
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Thus, the expansion of the Green's function is 

^''^^-^^=ikf2ji(kr^)h^\kry) J2 Yrrr,(e',^')YUO,^). (7.26) 



47r \r — r' 



/=0 m=-l 



This is a particularly useful result, as we shall discover, since it easily allows us to 
express the general solution of the inhomogeneous wave equation as a multipole 
expansion. 

It is well known in quantum mechanics that Eq. (7.6b) can be written in the 
form 

L^Yim = l{l + l)Yim. (7.27) 
The differential operator is given by 

= V + V + (7.28) 

where 

Z = -ir'AV (7.29) 
is times the orbital angular momentum operator of wave mechanics. 

The components of L can be conveniently written in the combinations 

L+ = L^ + iL^ = e^^(^^ + icot^^j, (7.30a) 

L_ = L,-iLy = e-'^ (^-^ + icote-^y (7.30b) 

L. = -iA. (7.30c) 

We note that L operates only on angular variables and is independent of r. From 
the definition (7.29) it is evident that 

r-L = (7.31) 

holds as an operator equation. It is easily demonstrated from Eqs. (7.30) that 

Id d 15^ 
= — -. — sin^— o— tt^- (7.32) 
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The following results are well known in quantum mechanics: 

L+Yim = ^/{l-m)(l + m + l)Yl,m+l, (7.33a) 



In addition, 



where 



L-Yim = ^/{l^m){l-m + l)Yi,m-i, (7.33b) 

LzYim = mYim- (7.33c) 

L'^L = LL^, (7.34a) 

L^L = iL, (7.34b) 

L^-V^ = V^Lj, (7.34c) 



V2 = ^— r^- ^. (7.35) 



J.2 Qj. Qj. J.2 

7.2 Multipole expansion of the vector wave equation 

Maxwell's equations in free space reduce to 



W-E = 


0, 


(7.36a) 


VcB = 


0, 


(7.36b) 


V AE = 


ikcB, 


(7.36c) 


V AcB = 


-ikE, 


(7.36d) 



assuming an e time dependence of all field quantities. Here, k = cu/c. Elimi- 
nating E between Eqs. (7.36c) and (7.36d), we obtain the following equations for 
B: 

(V^ + A;^)^ = 0, (7.37a) 

V B = 0, (7.37b) 

with E given by 

E= yV A cB. (7.38) 

k 
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Alternatively, B can be eliminated to give 

{V'^ + k'^)E = 0, (7.39a) 
V-^; = 0, (7.39b) 

with B given by 

cB = -yV AE. (7.40) 
k 

It is clear that each Cartesian component of B and E satisfies the Helmholtz 
wave equation (7.3). Hence, these components can be written in a general expan- 
sion of the form 

V(r) = Y, h^\kr) + 42 h^'Hkr)] y,„(«, v), (7.41) 

l,m 

where ip stands for any Cartesian component of E or cB. Note, however, that 
the three Cartesian components of E or B are not entirely independent, since 
they must also satisfy the constraints V--E = and V--B = 0. Let us examine 
how these constraints can be satisfied with the minimum labour. 

Consider the scalar r -A, where -4. is a well behaved vector field. It is easily 
verified that 

V'^{r-A) = r-(V^A) + 2VA. (7.42) 

It follows from Eqs. (7.37) and (7.39) that the scalars r- E and r- B both satisfy 
the Helmholtz wave equation: 

{V'^ + k'^){r-E) = 0, (7.43a) 
(V2 + A;2)(r-S) = 0. (7.43b) 

Thus, the general solutions for r-E and r cB can be written in the form (7.41). 

Let us define a magnetic multipole field of order (/,m) by the conditions 

^_^^(M) ^ ^l^g^^k^)YUe,^), (7.44a) 
vE\^^ = 0, (7.44b) 
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where 

gi {kr) = (kr) + /if ^ {kr) . (7.45) 

The presence of the factor /(/ + l)/k is for later convenience. Equation (7.40) 
yields 

kr-cB = -ir-(V A^;) = -i(r AV)-E = L E, (7.46) 

where L is given by Eq. (7.29). With r-B given by Eq. (7.44a), the electric field 
associated with a magnetic multipole must satisfy 

L-Ei^^ (r, e,ip) = l{l + 1) gi{kr) Yim{e, ^) (7.47) 

and r-Eil^^ = 0. Note that the operator L acts only on the angular variables 

{9,ip). This means that the radial dependence of e\^^ must be given by giikr). 
Note also, from Eqs. (7.33), that the operator L acting on Y/^ transforms the 
m value but does not change the I value. It is easily seen from Eqs. (7.27) and 
(7.31) that the solution to Eqs. (7.44b) and (7.47) can be written in the form 

e\^^ =gi{kr)LYim(eM- (7-48) 

Thus, the angular dependence of e'^^^ consists of some linear combination of 
Yi^rn-i, ^m, and Yi^rn+i- Equation (7.48), together with 



<=<' = -^VA<), (7.49) 



specifies the electromagnetic fields of a magnetic multipole of order {l,m). Note 
from Eq. (7.31) that the electric field given by Eq. (7.48) is transverse to the 
radius vector. Thus, magnetic multipole fields are sometimes termed transverse 
electric (TE) multipole fields. 

The fields of an electric or transverse magnetic (TM) multipole of order (/, m) 
are specified by the conditions 



r-E^^^ = --^^fi{kr)Yim(0,^), (7.50a) 



r-B^^^ = 0. (7.50b) 
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It follows that the fields of an electric multipole are given by 

c^lm = fi{kr)LYim{e,^), 

Im k 

The radial function fi{kr) is given by an expression like (7.45). 



(7.51a) 
(7.51b) 



The two sets of multipole fields (7.48), (7.49), and (7.51), form a complete set 
of vector solutions to Maxwell's equations in free space. Since the vector spherical 
harmonic LYim plays an important role in multipole fields, it is convenient to 
introduce the normalized form 



1 



(7.52) 



It can be demonstrated that the vector spherical harmonics possess the orthogo- 
nality properties 



j x;:,^,-{rAXim)dn = 0, 



mm' 5 



(7.53a) 
(7.53b) 



for all /, m, and m'. 



By combining the two types of fields we can write the general solution to 
Maxwell's equations in free space in the form 



l,m 



l,m 



ttE {I, m) fi (kr) Xim ~ -j^^M {I, m)V Agi (kr)Xi 



m 



(7.54a) 



- ttE {I, m) V A // {kr)Xim + Qm {I, m) gi (kr) Xi 



m 



(7.54b) 
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where the coefficients m) and aM{l, rn) specify the amounts of electric (/, m) 
and magnetic {l,m) multipole fields. The radial functions fi{kr) and gi{kr) are 
of the form (7.45). The coefficients «£■(/, m) and aM{li'rri), as well as the relative 
proportions in (7.45), are determined by the sources and the boundary condi- 
tions. 



Equations (7.54) yield 
rcB = 



T^ciM (I, m) gi (kr) L Xim 



k 

l,m 



(kr)^l(l + l)Yim, (7.55) 



and 



r E = 



k 

l,m 



- ^ X] m) fi {kr) L Xi 

l,m 



m 



= -lY.aE{l,m)fi{kr)^/l{[+l)Yirn, (7.56) 

l,m 

where use has been made of Eqs. (7.27), (7.29), and (7.31). It follows from the 
well known orthogonality property of the spherical harmonics that 

k f 

aM(l:m)gi{kr) = -j=== J Yi*^r-cB dfl, (7.57a) 

aE{l,m)fi{kr) = —-^=±== J Y^l^r-E dQ. (7.57b) 

Thus, knowledge of r-B and r-E at two different radii in a source free region 
permits a complete specification of the fields, including the relative proportions 
of /ip^ and /ip^ in // and gi. 
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7.3 Properties of multipole fields 



Let us examine some of the properties of the multipole fields (7.48), (7.49), and 
(7.51). Consider, first of all, the so-called near zone, for which kr <C 1. In this 
region fi{kr) is proportional to yi{kr), given by the asymptotic expansion (7.14b), 
unless its coefficient vanishes identically. Excluding this possibility, the limiting 
behaviour of the magnetic field for an electric {l,m) multipole is 

«<'^-yi^. (7-58) 

where the proportionality coefiicient is chosen for later convenience. To find the 
electric field we must take the curl of the right-hand side. The following operator 
identity is useful 

iVAZ = rV^- V^l + r^j . (7.59) 
The electric field (7.51b) is 

<)^^VAzf^V (7.60) 



Since Yim/r^'^^ is a solution of Laplace's equation, the first term in (7.59) vanishes. 
Consequently, the electric field at close distances for an electric {l,m) multipole 
is 

(7.61) 

This, of course, is an electrostatic multipole field. Such a field is obtained in a 
more straightforward manner by observing that E — )► — V(/), where = 0, in 
the near zone. Solving Laplace's equation by separation of variables in spherical 
polar coordinates, and demanding that be well behaved as |r| — >• oo, yields 

^(^.«.V') = E^W^- (7-62) 

l,m 

Note that the magnetic field (7.58) (normalized with respect to c~^) is smaller 
than the electric field (7.61) by a factor of order kr. Thus, in the near zone 
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the magnetic field associated with an electric multipole is always much smaller 
than the corresponding electric field. For magnetic multipole fields it is evident 
from Eqs. (7.48), (7.49), and (7.51) that the roles of E and B are interchanged 
according to the transformation 

e(^) ^ -cB^^\ (7.63a) 



In the so-called far zone or radiation zone, for which fcr 1, the multipole 
fields depend on the boundary conditions imposed at infinity. For definiteness, let 
us consider the case of outgoing waves at infinity, which is appropriate to radiation 
by a localized source. For this case, the radial function fi{kr) is proportional to 

the spherical Hankel function h^^\kr). From the asymptotic form (7.16), it is 
clear that in the radiation zone the magnetic field of an electric (/,m) multipole 
goes as 



cB^^'' ^ (-i) 



, i kr 



kr 



LY, 



Im • 



Using Eq. (7.51b), the electric field can be written 



(7.64) 



k^ 



ikr 



. i kr 



V 



/\LYim + 



V ALYi 



m 



(7.65) 



Neglecting terms which fall off faster than 1/r, the above expression reduces to 



i kr 



kr 



nALYim-^{rV'' -V)Yim 



(7.66) 



where use has been made of the identity (7.59), and n = r/r is a unit vector 
pointing in the radial direction. The second term in square brackets is smaller 
than the first term by a factor of order 1/kr, and can therefore be neglected in 
the limit /cr 1. Thus, we find that the electric field in the radiation zone is 
given by 

<' = MS A n, (7.67) 
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where Bi^^ is given by Eq. (7.64). These fields are typical radiation fields; i.e., 
they are transverse to the radius vector, mutually orthogonal, and fall off like 
1/r. For magnetic multipoles we merely make the transformation (7.63). 

Consider a linear superposition of electric {l,m) multipoles with different m 
values, but all possessing a common / value. It follows from Eqs. (7.54) that 



cBi = ^aE(l,m)Ximhl^\kr)e-'''\ (7.68a) 



El = yVAcBi. (7.68b) 
k 

For harmonically varying fields the time averaged energy density is given by 

u=j {E-E* + cB cB*). (7.69) 

In the radiation zone the two terms are equal. It follows that the energy density 
contained in a spherical shell between radii r and r -\- dr is 

1^ E aUl,m')aE{l,m) J X^^^-Ximdn, (7.70) 



m,m 



where the asymptotic form (7.16) of the spherical Hankel function has been used. 
Making use of the orthogonality relation (7.53a), we obtain 

''^ = ^Ei«^(''™)i'' ("1) 



dr 2A;2 



m 



which is clearly independent of the radius. For a general superposition of electric 
and magnetic multipoles the sum over m becomes a sum over / and m, and |a£;p 
becomes la^P + |£ImP- Thus, the total energy in a spherical shell in the radiation 
zone is an incoherent sum over all multipoles. 

The time averaged angular momentum density of harmonically varying elec- 
tromagnetic fields is given by 

m=^Re[rA{EA B*)]. (7.72) 
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For a superposition of electric multipoles the triple product can be expanded and 
the electric field (7.68b) substituted, to give 



(7.73) 



Thus, the angular momentum in a spherical shell lying between radii r and r-\-dr 
in the radiation zone is 



dM = 



egcdr 



Re^ a*E{l,m')aE{l,m) j {L-Xim')* XimdQ. (7.74) 

m,m' 



It follows from Eqs. (7.27) and (7.52) that 
dM eoc 



dr 2k^ 



Re^ a*E{l,m')aEil,m) j YC^, LYimd^. (7.7b) 

m,m' 



According to Eqs. (7.33), the Cartesian components of dM/dr can be written: 



dM^ 
dr 



dMy 
dr 



^Re^ [V(/-m)(/ + m+l)a|;(/,m+l) 

m 

+ \/ (/ + m)(l - m + 1) a^(/,m - 1) aE{l,m), 
^Im^ [V(^ - m)(/ + m + 1) a^il, m + 1) 



(7.76a) 



m 



■^y {I + m){l — m + 1) a*^{l, m — 1) aE{l,fn), (7.76b) 



dM, 
dr 



2P 



(7.76c) 



m 



Thus, for a general Ith. order electric multipole that consists of a superposition 
of different m values, only the z component of the angular momentum takes a 
relatively simple form. 
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7.4 Sources of multipole radiation 



Let us now examine the connection between multipole fields and their sources. 
Suppose that there exist localized distributions of electric change p{r,t), true 
current j{r, t), and magnetization ]VI{r, t). We assume that the time dependence 
can be analyzed into its Fourier components, and we therefore only consider 
harmonically varying sources, p{r) e~"^*, j{r) e~"^*, and M{r) e~"^*, where it is 
understood that we take the real parts of complex quantities. 

Maxwell's equations can be written 

V-^; = — , (7.77a) 

V B = 0, (7.77b) 
VAE-ikcB = 0, (7.77c) 
V AcB -\-ikE = /ioc(j + VAM), (7.77d) 

with the continuity equation 

ia;p = V-j. (7.78) 

It is convenient to deal with divergenceless fields. Thus, we use as the field 
variables, B and 

E' = E + ^j. (7.79) 

In the region outside the sources E' reduces to E. When expressed in terms of 
these fields. Maxwell's equations become 

V-^;' = 0, (7.80a) 
V B = 0, (7.80b) 

V AE' -ikcB = — VAj, (7.80c) 

V AcB -\-ikE' = /iocVAM. (7.80d) 

The curl equations can be combined to give two inhomogeneous Helmholtz wave 
equations: 

(V^ + k'^)cB = -i^ocV A (j + V A M), (7.81) 
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and 

(V^ + k'')E' = -i A; //oc V A (^M + ■ (7.82) 

These equations, together with V--B = 0, and V -E' = 0, and the curl equations 
giving E' in terms of B and vice versa, are the analogues to Eqs. (7.37) (7.40) 
when sources are present. 

Since the multipole coefficients in Eqs. (7.54) are determined according to 
Eqs. (7.57) from the scalars r -B and r ■ E' , it is sufficient to consider wave 
equations for these quantities, rather than the vector fields B and E'. From 
Eqs. (7.42), (7.81), (7.82), and the identity 

r-(V A A) = (r AV)-A = iL-A (7.83) 

for any vector field A, we obtain the inhomogeneous wave equations 

(V^ + A;2)r-cB = -i iiocL-{j + V A M), (7.84a) 

{V'' + k'')r-E' = kiiQcL-(M+^^\. (7.84b) 



Now the Green's function for the inhomogeneous Helmholtz equation (defined 
by Eq. (7.17) ), subject to the boundary condition of outgoing waves at infinity, 
is given by Eq. (7.18). It follows that Eqs. (7.84) can be inverted to give 

vcB{r) = ^ / L'-[j{r') + V AM{r')] dV, (7.85a) 



r-E'(r) = 



47r J 

k hqc 
47r 



\r — r\ 

/> „ i fe \r—r' I 



L' 



M{r') + 



V A j(r') 



(7.85b) 



In order to evaluate the multipole coefficients by means of Eqs. (7.57), we first 

(2) 

observe that the requirement of outgoing waves at infinity makes = in 

Eq. (7.45). Thus, we choose fi{kr) = gi{kr) = h^^\kr) in Eqs. (7.54) as the 
radial eigenfunctions of E and B in the source free region. Next, let us consider 



267 



the expansion (7.26) of the Green's function for the Helmholtz equation in terms 
of spherical harmonics. We assume that the point r lies outside some spherical 
shell which completely encloses the sources. It follows that r< = r' and r> = r in 
all of the integrations. Making use of the orthogonality property of the spherical 
harmonics, it follows from Eq. (7.26) that 

/ C\V-T'\ = ^ ^ ^^^^ ^^'^ ^^''^ (^-^^^ 

Finally, Eqs. (7.57), (7.85), and (7.86) yield 

a^(/,m) = -^^jMkr)Yr^L.(M+^^d'r, (7.87a) 

aM{l,m) = --^==Jji{kr)Yil,L-{j + VAM)d^r. (7.87b) 



The expressions (7.87) give the strengths of the various multipole fields outside 
the source in terms of integrals over the source densities j and M. They can be 
transformed into more useful forms by means of the following arguments. The 
results 

LA = iV-(rAA), (7.88a) 

L-(VAA) = iV'(r-A) -i-^^^^— ^ (7.88b) 

r or 

follow from the definition (7.29) of L, and simple vector identities. Substituting 
into Eq. (7.87a), we obtain 

aE(l,m} = --^0=Jj,(kr)Y,l,[V-(rAM) 



"' ' d^r, (7.89) 



kr dr 

where use has been made of Eq. (7.78). Use of Green's theorem on the second 
term replaces by —k^ (since we can neglect the surface terms, and ji{kr) 
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is a solution of the Helmholtz equation). A radial integration by part on the 
third term (again neglecting surface terms) casts the radial derivative over onto 
the spherical Bessel function. The result for the electric multipole coefficient is 



l^ock'^ f^^^ \ d[rji{kr)] •^ • /i, > 

cp — \-ik{r-j)ji{kr) 



-ikV-(rAM)ji(kr)]d'^r. (7.90) 

The analogous set of manipulations using Eq. (7.87b) leads to an expression for 
the magnetic multipole coefficient: 



2 



n \ /J-ock . 



V-(r A j) ji{kr) + V-M 



d[rji{kr)] 



dr 

- k^ (r-M) ji{kr)] d^r. (7.91) 

Both the above results are exact, and are valid for arbitrary wavelength and 
source size. 

In the limit in which the source dimensions are very small compared to a 
wavelength {i.e., kr <C 1) the expressions for the multipole coefficients can be 
considerably simplified. Using the asymptotic form (7.14a), and keeping only 
lowest powers in kr for terms involving p, j, and M, we obtain the approximate 
electric multipole coefficient 

where the multipole moments are 

Qim = jr'Yil^cpd^r, (7.93a) 

Q'lm = ^jr'YC^V-{rAM)d^r. (7.93b) 

The moment Qim has the same form as a conventional electrostatic multipole 
moment. The moment Q'l^ is an induced electric multipole moment due to the 
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magnetization. It is generally a factor kr smaller than the normal moment Qim- 
For the magnetic multipole coefficient aMiljfn) the corresponding long wave- 
length approximation is 

aM(l,m)cf^^ (^-±^y {Mim + MU, (7-94) 
where the magnetic multipole moments are 

Mim = -j^lr'Yil,V-irAj)d^r, (7.95a) 

M'lm = - j r'YC^V-Md^T (7.95b) 

Note that for a system with intrinsic magnetization the magnetic moments J^im 
and J^'i^ are generally of the same order of magnitude. 

Thus, in the long wavelength limit the electric multipole fields are determined 
by the charge density whereas the magnetic multipole fields are determined by 
the magnetic moment densities r A j/2 and M. 



7.5 Radiation from a linear cent re- fed antenna 

As an illustration of the use of a multipole expansion for a source whose di- 
mensions are comparable to a wavelength, consider the radiation from a linear 
centre-fed antenna. We assume that the antenna lies along the z-axis, and ex- 
tends from z = —d/2 to z = d/2. The current flowing along the antenna vanishes 
at the end points, and is an even function of z. Thus, we can write 

I{z,t)=I{\z\)e-'-\ (7.96) 

where I{d/2) = 0. Since the current flows radially, r Aj = 0. Furthermore, there 
is no intrinsic magnetization. Thus, according to Eq. (7.91), all of the magnetic 
multipole coefficients aM{h'm) vanish. In order to calculate the electric multipole 
coefficients aE{l,iTi), we need expressions for the charge and current densities. In 
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spherical polar coordinates the current density j can be written in the form 



j{r) = r 



27rr2 



[^(cos6'-l)-(5(cos6' + l)], 



(7.97) 



for r < d/2, where the delta functions cause the current to flow only upwards 
and downwards along the ^-axis. From the continuity equation (7.78), the charge 
density is given by 



p(r) = 



1 dl(r) 
iuj dr 



5(cose-l)-5{cose-\-l) 



(7.98) 



for r < d/2. 

These expressions for j and p can be substituted into Eq. (7.90) to give 



aE(l,'m) = 



/.d/2 



dr < kr ji{kr) I{r) — 



Idljr) djrjijkr)] 
dr 



} 



/ 



27rA//(/ + 1) Jo ' [ ' •"'^' ''^ ' k dr 

dn Yi*^ [5{cos e-1)- 5{cos 9 + 1)]. (7.99) 

The angular integral has the value 

J d^YC^ [5{cose- 1) - ^(cos^ + 1)] = 27r Ko(0) - Yi^i-n)] , (7.100) 

showing that only m = multipoles occur. This is hardly surprising given the 
cylindrical symmetry of the antenna. The m = spherical harmonics are even 
(odd) about ^ = 7r/2 for I even (odd). Hence, the only nonvanishing multipoles 
have / odd. So, the angular integral takes the value 



/ 



dnY{^[5{cosd - I) - 5{cose + I)] = v'47r(2/ + 1), (7.101) 
for / odd and m = 0. After some slight rearrangement, Eq. (7.99) can be written 



IIqc k 
27r 



47r(2/ + 1) 
/(/ + 1) 



1/2 



d 
dr 



rjiikr) — 



-\-rji {kr) 



dr"^ 



+ k^I 



dr. 



(7.102) 
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kd 




a£;(3,0)/a£;(l,0) 


a£(5,0)/a£(l,0) 




AVQtt (fiocI/4:7Td) 


4.95 X 10-2 


1.02 X 10-^ 


27r 


47r\/67r {hqcI / A'nd) 


0.325 


3.09 X 10-2 



Table 3: The first few electric multipole coefficients for a half- wave and a full- wave 
antenna 

In order to evaluate the integral (7.102) we need to specify the current I{z) 
along the antenna. In the absence of radiation, the sinusoidal time variation 
at frequency uj implies a sinusoidal space variation with wavenumber k = uj/c. 
However, the emission of radiation generally modifies the current distribution. 
The correct current I{z) can only be found be solving a complicated boundary 
value problem. For the sake of simplicity, we assume that I{z) is a known function; 
specifically, 

I(z) = I sin{kd/2 - k \z\), (7.103) 

for z < d/2, where / is the peak current. With a sinusoidal current the sec- 
ond term in curly brackets in Eq. (7.102) vanishes. The first term is a perfect 
differential. Consequently, Eqs. (7.102) and (7.103) yield 

(Y)%i(A'i/2), (7.104) 
for / odd. 

Let us consider the special cases of a half- wave antenna {kd = it; i.e., the 
length of the antenna is half a wavelength) and a full- wave antenna {kd = 27r). 
For these two values of kd the I = 1 coefficient is tabulated in Table 3, along with 
the relative values for / = 3, 5. It is clear from the table that the coefficients 
decrease rapidly in magnitude as / increases, and that higher I coefficients are 
more important the larger the source dimensions. However, even for a full-wave 
antenna it is generally adequate to retain only the / = 1 and / = 3 coefficients 
in order to calculate the angular distribution of the radiation. It is certainly 
adequate to keep only these two harmonics in order to calculate the total power 
radiated (which depends on the sum of the squares of the coefficients). 



aE{l,0) = 



7rd 



47t{21 + 1) 
/(/ + 1) 
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In the radiation zone the multipole fields (7.54) reduce to 



cB 



i(fer-wt) 



kr 



m 



l,m 



^aM(l,rn)n A Xim] , 



(7.105a) 

E ~ cBAn, (7.105b) 

where use has been made of the asymptotic form (7.16). The time-averaged power 
radiated per unit solid angle is given by 

dP Re(n-E AB*)r'^ 



2llQ 



(7.106) 



or 



dP 



dQ 2 fiQck' 



^(-i)'+^ [aE {1, m) Xim + ciM {I, m)nA Xim] 

l,m 



(7.107) 



Retaining only the / = 1 and / = 3 electric multipole coefficients, the angular 
distribution of the radiation from the antenna is given by 



dP \aE{lM' rv ^MMLrv 



(7.108) 



\/6as(l,0) 

where use has been made of Eq. (7.52). The various factors in the absolute square 
are 



1-^^3,0 1 



{LY^^oT-iLY^^o) = 



47r 
63 
167r 



sin2^(5cos2^- 1)2, 



sin2^(5cos2^-l). 



Stt 



(7.109a) 
(7.109b) 

(7.109c) 



With these angular factors, Eq. (7.108) becomes 

"7 0^(3,0) 



dP _ 3/ioc/2 3 sin2 6' 
— A 7, — 



dn 



Stt 



1 - 



8a£;(l,0) 



(5cos2 6'-l) 



(7.110) 
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where A equals 1 for a half- wave antenna and 7r^/4 for a full- wave antenna. The 
coefficient in front of {5cos'^9 - 1) is 0.0463 and 0.304 for the half- wave and 
full-wave antenna, respectively. It turns out that the radiation pattern from the 
two-term multipole expansion given above is almost indistinguishable from the 
exact result for the case of a half-wave antenna. For the case of a full-wave 
antenna the two-term expansion yields a radiation pattern which differs from the 
exact result by less than 5%. 

The total power radiated by the antenna is 

I odd 

where use has been made of Eq. (7.71). It is evident from Table 3 that a two-term 
multipole expansion gives an accurate expression for the radiated power for both 
a half- wave and a full- wave antenna. In fact, a one-term multipole expansion 
gives a fairly accurate result for the case of a half-wave antenna. 

It is clear from the above analysis that the multipole expansion converges 
rapidly when the source dimensions are of order the wavelength of the radiation. 
It is also clear that if the source dimensions are much less than the wavelength 
then the multipole expansion is likely to be completely dominated by the term 
corresponding to the lowest value of /. 

7.6 Spherical wave expansion of a vector plane wave 

In discussing the scattering or absorption of electromagnetic radiation by localized 
systems, it is useful to be able to express a plane electromagnetic wave as a 
superposition of spherical waves. 

Consider, first of all, the expansion of a scalar plane wave as a set of scalar 
spherical waves. This expansion is conveniently obtained from the expansion 
(7.26) for the Green's function of the scalar Helmholtz equation. Let us take the 
limit r' — ?► oo of this equation. We can make the substitution |r — r'| ~ r' — n-r 
on the left-hand-side, where n is a unit vector pointing in the direction of r' . On 
the right-hand side, r< = r and r> = r'. Furthermore, we can use the asymptotic 
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form (7.16) for h\^\kr). Thus, we obtain 



= '*^^rEH)''''.''(^'-)^™(«''V')y/™(«,v)- (7.112) 

l,m 

Canceling the factor e''^^ /r' on either side, and taking the complex conjugate, 
we get the following expansion for a scalar plane wave, 

oo +1 

e'fc- = 47r^i'i,(fcr) ^ Y,l,ie,^)Yim{e' ,^'), (7.113) 

/=0 m=-l 

where k is the wave vector with the spherical coordinates k, 9', cp'. The well 
known addition theorem for the spherical harmonics states that 

^'('°'^^=2lTl ^ yrm{0,^)Yim(0',^% (7.114) 

m=— 1 

where 7 is the angle subtended between the vectors r and r'. Consequently, 

cos 7 = cos^ cos^' + sin^ sin^' cos{cp — cp'). (7.115) 
It follows from Eqs. (7.113) and (7.114) that 



00 



^ife.r ^ (2l + l)ji(kr)Pi(cos^), (7.116) 



1=0 

or 

00 



e*'^-^ = ^iV47r(2/ + l)iz(fcr)11,o(7), (7-117) 



since 



no(S) = ^^P,(cosO). (7,118) 

Let us now make an equivalent expansion for a circularly polarized plane wave 
incident along the z-axis: 

E{r) = {x±iif)e'^', (7.119a) 
cB{r) = zAE = ^iE. (7.119b) 
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Since the plane wave is finite everywhere (including the origin), its multipole 
expansion (7.54) can only involve the well behaved radial eigenfunctions ji{kr). 
Thus, 



E = 



cB = 



E 

l,m 



E 

l,m 



a± {I, rn) ji (kr) Xim + ^ ^± m) V A ji {kr)Xi 



m 



(7.120a) 



— 1 

T 



a± {I, m) V A ji {kr)Xim + h± {I, m) ji {kr) Xim 



(7.120b) 



To determine the coefficients a±{l,m) and b±{l,i7i) we make use of a slight gen- 
eralization of the standard orthogonality properties (7.53) of the vector spherical 
harmonics: 

J[fi{r)Xi>m'T ■[9i{r)Xim]dn = f^gi6ii>5mm', (7.121a) 
Jlfi(r)Xi>m'T -[^ Agi(r)Xim]dn = 0. (7.121b) 

The first of these follows directly from Eq. (7.53a). The second follows from 
Eqs. (7.31), (7.53b), (7.59), and the identity 



„ r d i 

V = ^rAL. 

r or r^ 



(7.122) 



The coefficients a±{l,m) and b±{l,m) are obtained by taking the scalar prod- 
uct of Eqs. (7.120) with Xi^ and integrating over all solid angle, making use of 
the orthogonality relations (7.121). This yields 

a±{l,m)ji(kr) = J X^^-EdQ, (7.123a) 
b±(l,m)ji(kr) = J Xtm-cBdQ. (7.123b) 
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Substitution of Eqs. (7.52) and (7.120a) into Eq. (7.123a) gives 



(7.124) 



where the operators L± are defined in Eqs. (7.30). Making use of Eqs. (7.33), the 
above expression reduces to 

a^il,m)Mkr) = ^^J^ + J YC^,, e^'^Un. (7.125) 

If the expansion (7.117) is substituted for e^^^, and use is made of the orthogo- 
nality properties of the spherical harmonics, then we obtain the result 



(7.126) 



a±{l,m) = i' a/^tt (2/ + 1) 
It is clear from Eqs. (7.119b) and (7.123b) that 

b±{l,m) = ^ia±{l,m). (7.127) 
Thus, the general expansion of a circularly polarized plane wave takes the form 



oo 



E{r) = ^iV47r(2/ + l) 



1=1 



oo 



jiikr)Xi,±i ± - V Aji{kr)Xi,±i 



(7.128a) 



B{r) = ^iV47r(2/ + l) 



1=1 



— 1 



— V A ji{kr)Xi^±i ^ iji{kr)Xi^±i 



(7.128b) 



The expansion for a linearly polarized plane wave is easily obtained by taking the 
appropriate linear combination of the above two expansions. 



7.7 Mie scattering 



Consider a plane electromagnetic wave incident on a spherical obstacle. In gen- 
eral, the wave is scattered, to some extent, by the obstacle. Thus, far away from 
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the sphere the electromagnetic fields can be expressed as the sum of a plane wave 
and a set of outgoing spherical waves. There may be absorption by the obstacle, 
as well as scattering. In this case, the energy flow away from the obstacle is 
less than the total energy flow towards it: the diffierence represents the absorbed 
energy. 

The fields outside the sphere can be written as the sum of incident and scat- 
tered waves: 

E{r) = Einc + Esc, (7.129a) 
Bir) = Bi^^^Bsc, (7.129b) 

where E^^c and -Bine are given by (7.128). Since the scattered fields are outgoing 
waves at infinity, their expansions must be of the form 



Esr 



2 

1=1 



cB^r = 



-^iV47r (2/ + 1) (7.130a) 

a± (/) (kr) ± Ml V A (kr) 

-^iV47r (2/ + 1) (7.130b) 

i=\ 

V A hf\kT)X,^^^ ^ i/3±(0 hf\kT) 

The coefl&cients ol±{J) and ^±{p) are determined by the boundary conditions on 
the surface of the sphere. In general, it is necessary to sum over all m har- 
monics in the above expressions. However, for the restricted class of spherically 
symmetric scatterers only m = ±1 harmonics need be retained (since only these 
harmonics occur in the spherical wave expansion of the incident plane wave (see 
Eqs. (7.128) ), and a spherically symmetric scatterer does not couple different m 
harmonics). 

The angular distribution of the scattered power can be written in terms of the 
coefficients a(/) and /5(/) using the scattered electromagnetic fields evaluated on 
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the surface of a sphere of radius a surrounding the scatterer. In fact, it is easily 
demonstrated that 



dR 



sc 



a 



2/io 



a 



2/20 



Re [E,^-(nAB;^)]r=a, 



(7.131) 



where n is a radially directed outward normal. The differential scattering cross 
section is defined as the ratio of dP^ddQ to the incident flux Ij hqc. Hence, 



^ = -yRe[£;,e-(nAcB;j],=„. 



(7.132) 



We need to evaluate this expression using the electromagnetic fields specified in 
Eqs. (7.128), (7.129), and (7.130). The following identity, which can be estab- 
lished with the aid of Eqs. (7.29), (7.52), and (7.59), is helpful in this regard: 



r r dr 

For instance, using this result we can write n A cBsc in the form 

^ oo 

AcBsc = -^iV^TT (2/ + 1) 

a±{l) ld[rhl^\kr)] 



(7.133) 



n 



1=1 



(7.134) 



k 



dr 



It can be demonstrated, after considerable algebra, that 
d<T. 



'sc 



TV 



dQ 2A;2 



^ V'2Hn:[a±(0 ± i/3±(/) n A 



(7.135) 



In obtaining this formula, use has been made of the standard result 

Ji'\z) :n— jM = 72' 

At 



dz 



dz 



(7.136) 
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where fi{z) = h^^\z). The total scattering cross section is obtained by inte- 
grating Eq. (7.135) over all solid angle, making use of the following orthogonality 
relations for the vector spherical harmonics (see Eqs. (7.53) ): 



/ 



/ 



^l'm''^lmdfi — Sii'dmm', (7.137a) 

j Xl^,\nAXim)dn = 0, (7.137b) 
{nAX;;,^,)-{nAXim)dn = SwSmm'- (7.137c) 



Thus, 



According to Eqs. (7.135) and (7.138), the total scattering cross section is inde- 
pendent of the polarization of the incident radiation [i.e., it is the same for both 
the ± signs). However, the differential scattering cross section in any particular 
direction is, in general, different for different circular polarizations of the incident 
radiation. This implies that if the incident radiation is linearly polarized then the 
scattered radiation is elliptically polarized. Furthermore, if the incident radiation 
is unpolarized then the scattered radiation exhibits partial polarization, with the 
degree of polarization depending on the angle of observation. 

The total power absorbed by the sphere is given by 

.2 



Pahs — 



— Re / [n-EAB*]r=a dQ 

-Re f [E-(nAB*)]r=adn. (7.139) 
J 



A similar calculation to that outlined above yields the following expression for 
the absorption cross section, 

^-bs = ^ + 1) [2 - l«±(0 + - + IP] ■ (7.140) 
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The total or extinction cross section is the sum of osc and (Tabs'- 



(7.141) 



Not surprisingly, the above expressions for the cross sections closely resemble 
those obtained in quantum mechanics from partial wave expansions. 

Let us now consider the boundary conditions at the surface of the sphere 
(whose radius is a, say) . For the sake of simplicity, let us suppose that the sphere 
is a perfect conductor. In this case, the appropriate boundary condition is that 
the tangential electric field is zero at r = a. According to Eqs. (7.128), (7.129), 
and (7.133), the tangential electric field is given by 



E 



tan 



5]iVM2/ + i)| 



ji + 



a±(l) 



X 



i,±i 



^ 1 d 

X dx 



n A Xi^±x 



(7.142) 



where x = ka, and all of the spherical Bessel functions have the argument x. 
Thus, the boundary condition yields 



a±(/) + l = 



/5±(/) + l = 



{^hf\x)y 
{xh\'\x)y 



(7.143a) 



(7.143b) 



x=ka 



where ' denotes d/dx. Note that a±{l) + 1 and [3±{l) + 1 are both numbers of 
modulus unity. This implies, from Eq. (7.140), that there is no absorption for the 
case of a perfectly conducting sphere (in general, there is some absorption if the 
sphere has a finite conductivity). We can write ol±{1) and [3±{l) in the form 



a±(/) = e^'^'-l. 



(7.144a) 
(7.144b) 
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where the phase angles Si and d'l are called scattering phase shifts. It follows from 
Eqs. (7.143) that 



tan 6i 



tan 6', = 



ji(ka) 
yi{kay 

(xji(x)y 

{xyi{x) y 



(7.145a) 
(7.145b) 



x=ka 



Let us specialize to the limit fca <C 1 , in which the wavelength of the radiation 
is much greater than the radius of the sphere. The asymptotic expansions (7.14) 
yield 

a±{l) ~ 



(2/ + l)[(2/-l)!!]2' 

P±(l) ^ -^^^a±(0, (7.146a) 

for / > 1. It is clear that the scattering coefficients a±{l) and (3±{l) become small 
very rapidly as / increases. In the very long wavelength limit only the / = 1 
coefficients need be retained. It is easily seen that 

«±(l) = -^--yM^ (7.147) 
In this limit, the differential scattering cross section (7.135) reduces to 

"^^f a\kaf |JCi,±i T 2i « A JCi,±i ^ • (7-148) 

It can be demonstrated that 

|n A = = -|- (1 + cos^ 9), (7.149) 

and 

Stt 

[±i(n AX*±i)-Xi,±i] = -— cos^. (7.150) 
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Thus, in long wavelength limit the differential scattering cross section limits to 



dn 



a^{kaY -(l + cos^6l)-cos^ . 



(7.151) 



The scattering is predominately backwards, and is independent of the state of 
polarization of the incident radiation. The total scattering cross section is given 



This well known result was first obtained by Mie and Debye. Note that the 
cross section scales as the inverse fourth power of the wavelength of the incident 
radiation. This scaling is generic to all scatterers whose dimensions are much 
smaller than the wavelength. In fact, it was first derived by Rayleigh using 
dimensional analysis. 



by 




(7.152) 
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